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The equations for phase transitions temperatures, order parameters and critical concentrations 
of components have been derived for mixed ferroelectrics. The electric dipoles randomly distributed 
over the system were considered as a random field sources. We derive a random field distribution 
function for different orientations of the electric dipoles with nonlinear and spatial correlation ef- 
fects included. The concentrational dependence of transition temperature for Ai-iB^! mixed system 
has been calculated for A and B being antiferroelectric and ferroelectric materials as well as for A 
and B being ferroelectric and paraelectric. The numerical calculations have been carried out for 
PbZri-i^Tii^Oa and BaZr^jTii-ajOs. The obtained results give a fairly good description of exper- 
imentally observed phase diagrams of these mixed systems. We discuss the physical reasons of 
strongly different behaviour of the systems related (in particular) to the relaxor properties appear- 
ance in BaZr^jTii-ajOs at a; > 0.27 and peculiar role of lead ions. We predict a transformation 
of any mixed system with ferroelectric (antiferroelectric) and paraelectric component into relaxor 
material at some large enough concentration of paraelectric component. 



I. INTRODUCTION 



There has been considerable amount of attention of sci- 
entists to mixed systems of different types like ferroelec- 
tric and antiferroelectric, ferroelectric and paraelectric 
etc due to their unusual physical properties, utilizable in 
many applications. A classic example of such a system 
is PbZri_2;Tia;03 (PZT) with broad range of its tech- 
nical usage [y, including novel branches of electronics 
Q . Most of the important peculiarities occur in the spe- 
cial regions of concentrations in the vicinity of different 
phase boundaries on phase diagram. In particular, for 
PZT the most interesting range is the morphotropic one 
with coexistence of ferroelectric phases with two differ- 
ent symmetries and antiferroelectric-ferroelectric phase 
boundary pjlO[ . The characteristic features of phase di- 
agram and properties anomalies depend essentially on 
the type of solid solution components (ferroelectric, an- 
tiferroelectric or paraelectric). It is especially interest- 
ing from this point of view to compare phase diagrams 
of PZT (solid solution of antiferroelectric PbZrOs and 
ferroelectric PbTiOs) and BaZr^Tii_^03 (BZT, which 
is solid solution of ferroelectric BaTiOs and paraelectric 
BaZrOs |14|. The latter clearly demonstrates a peculiar 
role of lead ions in phase transitions of perovskite struc- 
ture materials (see e.g. [M). Due to different states (e.g. 
ferroelectric, paraelectric etc) of PZT and BZT compo- 
nents, their phase diagrams are strongly different (see 
e.g. [OJ3]). Namely, several phases with long range or- 
der exist in PZT while BZT exhibits relaxor behaviour 



beginning at a; = 0.27. To find out the physical mech- 
anisms "responsible" for actual phase diagram forma- 
tion in mixed system, the theoretical calculations seem 
to be extremely desirable. Recently, the random field 
model for calculation of physical properties of mixed sys- 
tem had been proposed (l7|. In the present work this 
model is expanded by taking into account the nonlinear 
and correlation effects as well as different orientations 
of the dipoles. This made it possible to calculate the 
concentration dependence of transition temperature, de- 
termine the symmetry of order parameters (and obtain 
the change of mixed system symmetry) . This model also 
permits to obtain the range of parameters for existence of 
morphotropic region and appearance of glassy state. We 
apply this essentially improved model to calculations of 
phase diagrams of PZT and BZT. The theory describes 
adequately the observed phase diagram in these materi- 
als. 



II. ORDER PARAMETERS 

In the case of solid solutions of antiferroelectric and 
ferroelectric materials we have to consider three order 
parameters. Namely, in two- sublattice model for antifer- 
roelectric component there are ferroelectric L2F and anti- 
ferroelectric L2A order parameters that describe, respec- 
tively, the homogeneous and inhomogeneous displace- 
ments of ions. The third order parameter Lip is related 
to another component of solid solution. Keeping in mind 
the key role of electric dipoles in the system and their 



ability to be oriented by electric field E we write the or- 
der parameters in the form inherent for Ising model |l8[ : 
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Here d\ and d^ are effective dipole moments related to 
the first (a; = 1) and the second [x = 0) components; Tip 
is ferroelectric phase transition temperature of the first 
component and T2F, T2A are the transition temperatures 
of ferroelectric and antiferroelectric phases of the second 
component of solid solution. Superscript " mf " in L shows 
that Eqs.(l) are written in a mean field approximation. 
In the mixed system the positions of the dipoles of both 
components become random so that they can be consid- 
ered as a random field sources. In this case order pa- 
rameters given by Eqs. (1) have to be averaged with the 
random field distribution function f{E, Lip, L2f)- Note, 
that in general case electric field E is the sum of external 
and internal field, but hereafter we suppose the absence 
of external field. Allowing for nonlinear and spatial corre- 
lation effects contribution into the distribution function, 
H order parameters for the mixed system can be written 
as: 
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Here ipi{E) is the nonlinear function of the field, which 
has the form of infinite series in odd powers of E in the 
lattices with a center of inversion in paraelectric phase. 
We keep in (Gq) only first nonvanishing nonlinear term 



with the coefficient as. Also, f{E,Lip,L2p) is a lin- 
ear random field distribution function. It is seen, that 
in the case when both components are ferroelectrics, i.e. 
L2A — 0, Eq.(2b) has the same form as Eq.(p3) as it has 
to be expected, but with different parameters. There- 
fore Eqs. (2) describe a quite general case. They give 
the order parameters dependence on the characteristics 
of mixed system components and their molar fractions 
via the random field distribution function. We would 
like to emphasize that the parameters L are the fraction 
(0 < -^ < 1) of coherently oriented dipoles of the mixed 
system components. Polarization P as the actual order 
parameter of mixed system can be expressed via Li as: 
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where oi, 02 and x, 1—x are respectively lattice constants 
and molar fractions of the first and the second compo- 
nents in the solid solution with chemical formula Ai-^jB^:. 
Note, that L2A contributes neither to mixed system po- 
larization nor to the random field distribution function as 
it has to be expected. The orientation of mixed system 
polarization is related to the vector sum of the dipoles in 
accordance with Eq.(g). To derive this sum coefficients, 
one has to calculate Li,2P on the base of Eqs. (2) which 
depend on the form of function f{E, Lip, L2f). 



III. THE RANDOM FIELD DISTRIBUTION 
FUNCTION 



As we have seen above, to calculate the nonlinear dis- 
tribution function of random fields, it is sufficient to cal- 
culate the linear distribution function (i.e. that without 
nonlinear and spatial correlation effects) . For latter pur- 
pose all kind of dipoles can be considered as independent 
sources of random field. Therefore, the distribution func- 
tion of the mixed system is a convolution of the two types 
of dipoles distribution functions B. Gaussian approxi- 
mation for these functions leads to the following expres- 
sion for mixed system distribution function: 
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Here Eg is a mean field, A_Ei and /S.E2 are the 
halfwidths of the distribution functions induced by d\ 
and ^2 dipoles respectively. They can be written in the 



form 0: 



where £1 and rd, £2 and rrj2 are dielectric permittivities 
and correlation radii of the solid solution components. 

One can see from Eqs. (H), (p[), (ra) that mixed system 
distribution function depends on the components concen- 



trations, order parameters Lip L2F, transition tempera- 
tures Tip T2F, their dipole moments and other physical 
parameters. 



IV. PHASE DIAGRAM. GENERAL EQUATIONS. 

The phase diagram has to reflect the variation of the 
transition temperatures of different phases with the in- 



creasing of the fraction of the components as wefl as the 
phases symmetry changing. All this information can be 
obtained by solving the equations (2) with respect to Eqs. 
(|), (|). Substitution of Eq.(|) into (2) leads to six-fold 
integrals. They can be simplified due to the dependence 
of all the integrands on the scalar product (Ee.i), i.e. only 
on the parallel to e.^ field component. This permits to in- 
tegrate out two other field components and to perform 
the integration over cfip. This finally yields: 
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where i?oi — (Eoe^) [i = 1, 2), so that E'oi depends on the 
angle between the directions of the two types of dipoles, 
because (6162) = cos(ei,e2) = cos9. The parameters A^ 
also depend on this angle, namely: 



Ai ^ x{l^Eif + {I - x){l^E2C0se 
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The Eqs.(0) are the final form of the equations for the 
order parameters L2A, L2f, Lip dependence on the mo- 
lar fraction x and via them polarization of the system (see 
Eq.(||)). The transition temperature Tc to the ferroelec- 
tric phase and the transition temperature Ta to antifer- 
roelectric one in the mixed system can be derived from 
Eq.(7) in the limit of zeroth order parameters. One can 
obtain from Eq.(7) the following system of equations: 
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Here dimensionless variables are introduced: 
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One can see that the temperature Ta follows from Eq.(9c), while the temperature Tc can be derived from Eqs. (9a), 
(9b). The solution of these equations has the form: 
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The solution of Eqs (12) gives the dependences of Tc 
and Ta on molar fractions and material parameters of 
the mixed system components (see (11)). 



V. COMPARISON OF THE THEORY AND 
EXPERIMENT 

A. Phase diagram of PbZri-ajTi^^Os 

This solid solution components are antiferroelectric 
PbZrOs with transition temperature T2A = 503 K (the 
value T2F ~ T2A @) and ferroelectric FbTiOs with 
transition temperature from paraelectric phase to tetrag- 
onal ferroelectric phase Tip — 763 K. Both components 
have electric dipoles randomly distributed in the mixed 
system. In accordance with the components symmetry 
one can suppose that d\ \\ [001] and d^ \\ [HI] types of 
directions. Therefore parameters Xp ~ 1.516, A^ = 1 
and cos 6* — l/\/3. Other parameters were obtained from 
the fitting with observed phase diagram of PZT. We be- 
gin with the fitting of Ta from Eq.(12b). The numerical 
calculations of the integral I2 were performed at a2 = 0.3, 
(72 = 2.9. The dependence of Ta on molar fraction x is 
depicted by dashed line in Figure 1. 
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FIG. 1. Phase diagram of PbZri_i,Tia;03. Open squares 
are the experimental data |I3]. Calculated transition tem- 
peratures are the following: solid line represents the transi- 
tion from paraelectric phase (Pc) to ferroelectric phases (Fr, 
Ft), dashed line represents the transition from antiferroelec- 
tric phase (AF) to ferroelectric rhombohedral phase (Fr), dot- 
ted line represents the transition from (Fr) phase to ferroelec- 
tric tetragonal phase (Ft). The transition between Pc and Fr 
phases at small fractions x is depicted in inset. 

One can see that Ta decreases with x increase and 
Ta = K at a; = ccc = 0.093 (where xc is a critical frac- 



tion at which antiferroelectric phase disappears). It is 
seen that the theory gives reasonable fit to experimental 
data shown by open squares in Figure 1. The transition 
temperature Tc calculated on the base of Eq.(12a) in- 
creases with X increase (see solid line in Figure 1). Our 
theory is undoubtedly valid at small x (see inset to Figure 
1) . The fitting of Tc for all the range of x were performed 
by varying of p, 51, 52, Oi, 012- The best fit was achieved 
at p = 0.828, qi = 0.239, ga = 0.364, ai = 3.9, as = 
4.3. As it follows form Figure 1, the fitting is also good 
for x > 0.6, i.e. in the region enriched by titanium. For 
intermediate molar fractions 0.1 < a: < 0.6 the accuracy 
of the fitting is not so good as it is in the other regions. 
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FIG. 2. Phase diagram BaZriTii-iOs. Open circles, 
squares and triangles are the experimental data from S. 
Calculated transition temperatures are the following: dashed 
line represents the transition between tetragonal (Ft) and or- 
thorhombic (Fo) ferroelectric phases, dotted line represents 
the transition from (Fo) phase to ferroelectric rhombohedral 
phase (Fr), solid line represents the transition from paraelec- 
tric cubic phase (Pc) to (Ft, Fr) phases. 



To our mind this deviation is related to the existence 
(in this interval of x)of additional order parameters orig- 
inated from an improper ferroelectric phase transition. 
This transition occurs between high and low temperature 
rhombohedral phases. In the high temperature phase the 
spontaneous tilting of the oxygen octahedra causes the 
phase transition and contributes to the spontaneous po- 
larization uM- The symmetry of different ferroelectric 
phases and morphotropic region with both symmetries 
coexistence (see dotted line in Figure 1) was calculated 
on the base of Eqs. (3), (7), (8) with the same set of 
parameters as that used for the transition temperatures 
fitting. It is clear from Figure 1 that the morphotropic 
region lies between x — 0.453 at T = K and x = 0.463 
at T = 611 K. Therefore the calculations describe ade- 



quately the observed phase diagram Q represented by 
open squares in the Figure 1. Note that the best fit was 
obtained in the assumption that the halfwidths A^ of the 
distribution functions for each component are determined 
by the other component of sohd solution. This confirms 
the supposition that the sources of the random field of 
the first component destroy the long- range order of the 
second component and vise versa. 



B. Phase diagram of BaZr^^Tii-iiOs 

The main component of this solid solution is the fer- 
roelectric material BaTiOa which is known to have three 
phase transitions at T^]^ = 400 K, T^'^J = 305 K, 
T^p = 184 K, to ferroelectric phases with tetragonal, or- 
thorhombic and rhombohedral symmetries respectively. 
BaZrOs is paraelectric at all temperatures [O. Phase 
diagram of BZT differs strongly from that of PZT. To 
describe it we proposed the following model. 

We assume that in mixed BZT system zirconium ions 
can be shifted so that they supply the random electric 
dipoles to BaTiOs component. These dipoles are the 
main sources of the random field. This field distribution 
function halfwidth has to be larger than the mean field 
produced by zirconium dipoles because BaZrOs is para- 
electric. Taking all these arguments into consideration 
we describe the phase diagram of BZT system with the 
help of Eqs.(3), (7), (8) with the following set of param- 
eters: 

TiF ^ 250 K, qi = 0.6, 92 =0, ai= 0; 
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= 2, p(2) = 1.63; (13) 



Note, that is would be more accurate to write Ti,„/ in- 
stead of Tip because there is no actual ferroelectric phase 
transition in BaZrOs. 

Allowing for known symmetry of three BaTiOa ferro- 
electric phases, one can obtain cos(6''^^^) = 1, cos(6'*^^)) = 
I/V2, cos(6'(3)) = l/\/3. The results of the calculation 
are shown in Figure 2. One can see that the theory 
describes the observed phase diagram (see |^ and ref- 
erences therein) quite well. Note, that the accuracy of 
fitting of experimental points by dashed and dotted lines 
at a; > 0.12 is about 10%. 

In Figure 3 we represented BZT order parameters cal- 
culated with the parameters (13) at T = K. It is easy 
to check that in this limit ii = L2. One can see that 
at X > 0.3 the fraction of coherently oriented dipoles L 
is less than 0.9 which corresponds to mixed ferroglass 
phase H, the critical fraction for the dipole glass ap- 
pearance [L = 0) being about Xc ~ 0.82. The relaxor 
behaviour (e.g. Vogel-Fulcher law in dynamic permit- 
tivity) was observed recently at x > 0.27 (see p and 
references therein). 



The value of Xc for the transition to dipole glass state 
is a prediction of the theory. Unfortunately, the avail- 
able experimental data cover the range up to a; = 0.5. 
The clarification of dipole glass state existence may be 
performed in solid solution of BaTiOs and BaZrOa only 
under the condition of the components solubility in all 
the concentration range. 
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FIG. 3. The dependence of BaZra;Tii_^03 order parame- 
ters at T = K versus the molar fraction x. 



VI. DISCUSSION 



In the proposed model two components of the mixed 
system were considered as the materials with electric 
dipoles di and d2. The dipoles tend to order the sys- 
tem along their directions so that the competition be- 
tween different ordering directions is the main mechanism 
defining phase diagram of mixed system. Without this 
competition (when there is only one type of dipoles) it 
can be ferroelectric (PbTiOa, BaTiOa), antiferroelectric 
(PbZrOs) or paraelectric (BaZrOa) phases. The differ- 
ence between properties of the materials with Pb or Ba 
ions is related to peculiar role of lead ions in the phase 
transitions. In perovskite structure ABO3 all A ions, but 
Pb, give almost no contribution to a lattice polarization. 
In contrast to this Pb ions contribution (e.g. in PbTiOa) 
is the main one (see e.g. ||l^ and ref. therein). Oppo- 
site displacements of lead ions in PbZrOa are known to 
be the characteristic feature of antiferroelectric phase in 
PbZrOa. In our model d^ represents the resultant dipole 
moment of a lattice unit cell. It is a vector sum of Pb 
and Ti ions displacements in the case of PbTiOa or Ti 
ions in BaTiOa, all the displacements being considered 
relatively oxygen cage. Therefore the extraction of lead 
contribution can be made only on the base of independent 
microscopic calculations or measurements of d\/d2 ratio 
which is fitting parameter in our model. The estimation 
of this ratio with the help of known displacements (see 
e.g. fll) of Ti ions in PbTiOa and Zr ions in PbZrOa 



gives p = dydl « 0.8 while our fit gives p = 0.828 (see 
section 5.1). Since nothing is known about Zr displace- 
ments in BaZrOs it is resonable to calculate, e.g., the 
ratios of fci = pi/pa and k2 = P2/P3, which are indepen- 
dent on dipole moment of Zr ions. Under supposition 
that fci and fc2 can be estimated via the ratio of polariza- 
tions in ferroelectric phases of BaTiOa (see e.g. [Q), one 
obtains ki — 1.8, fc2 = 1.4 while in our model fci = 1.6, 
^2 = 1.3 (see (13)). The values oi p in (13) show that Zr 
displacement in BaZrOs is two times smaller, than that of 
Ti in tetragonal phase of BaTiOa. Therefore our model 
leads to reasonable values of the ions dipole moments 
both in PZT and BZT. It is actually possible to estimate 
fitting parameters with the help of Eqs.(6) provided that 
correlation radii and dipole moments are known. The 
parameters related to nonlinearity coefficients play an 
important role in description of peculiar form of BZT 
phase diagram because they are "responsible" for Tc{x) 
maxima. Although the relation between ckj is qualita- 
tively the same as that between nonlinearity constants a 



in different phases of BaTiOa (namely a^' < a^ < a 
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), quantitative es- 
timation of ttj''' seems to be complex problem and the 
independent measurements of these parameters are de- 
sirable. Therefore these parameters are the actual fitting 
ones while other parameters could be estimated a priori. 
Another interesting feature of BZT phase diagram is ap- 
pearance of relaxor behaviour at x > 0.27 while in PZT 
nothing of this kind is known. This may be related to 
the fact that in PbTiOa (due to lead ions displacements) 
spontaneous polarization is more than 2 times larger than 
that in BaTiOa |lq| . As a result the random field induced 
by Zr ions appeared unable to destroy strong ferroelectric 
order in PbTiOa whereas in BaTiOa this field can destroy 
weaker ferroelectric long range order at large enough Zr 
ions concentration. 

More generally we can assert that for a mixed system 
containing ferroelectric and paraelectric components, it 
has to be the concentration range where the system trans- 
forms into relaxor. This statement follows from the fact 
that paraelectric component contribution to mean field 
Eq is rather small while it completely defines the dis- 
tribution function half-width (92 = 0, see (13)). The 
contribution of the ferroelectric component to Eq de- 
creases (see Eq. (5) with subscripts 1 and 2 correspond- 
ing to paraelectric and ferroelectric component respec- 
tively and T2A — 0, L2A = in Eqs. (7)). So, the in- 
crease of paraelectric component concentration leads to 
mean field decrease and to increase of distribution func- 
tion half-width A (see Eqs. (6), (8)). This must result 
into Eo/A{x) decrease. In supposition that the state of a 
system (paraelectric (PE), ferroelectric (FE), dipole glass 
(DG), mixed ferroglass (EG) where FE long range order 
coexists with DG short range order) strongly depends on 
ratio Eq/A (see Figure 4), one can conclude that as x 



increases, the system passes from FE to EG and than 
to DG state at some low temperature region (see arrows 
in Figure 4). Both EG and DG states are known to be 
characteristic feature of relaxor materials [Q. The re- 
laxor systems exhibit nonergodic behaviour just in these 
phases. Vogel-Fulcher law which describes temperature 
dependence of dynamic dielectric susceptibility of relax- 
ors can be related to the distribution of random electric 
fields in the mixed system m. The calculations of con- 
centrational dependence of (BaTi03)i_a;(BaZr03)a; or- 
der parameter made it possible to obtain the concen- 
tration (x w 0.3, see Figure 3) at which mixed system 
transforms into relaxor. It follows from Figure 3 that 
BZT is in EG phase at 0.3 < a; < 0.8 and could be 
in DG state at x > 0.8, the latter being dependent on 
existence of this x region in BZT mixed system. The 
existence of DG state and relaxor properties in another 
mixed system (BaTi03)i-x(SrTi03):r (BST) at x > 0.9 
|l5| confirms the generality of the statement about trans- 
formation of mixed ferroelectric-paraelectric system into 
relaxor for some concentration range. This transforma- 
tion seems to exist also in the mixed system consisting 
of antiferroelectric (component 2) and paraelectric (com- 
ponent 1). Really, from Eqs. (5), one can expect the de- 
crease of Eq/A so the system at x increase passes from 
FE to EG and DG states (see the arrow in Figure 4). 
However, contrary to the case of ferroelectric-paraelectric 
mixed system, where T2F is actual temperature of ferro- 
electric phase transition, in the above considered case 
T2F is characteristic of "imaginary" ferroelectric phase 
following from two sublattice model of antiferroelectrics. 
Its value is close (although little lower) to T2A ll^ • 
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FIG. 4. Scheme of the disordered system phase diagram B 



VII. CONCLUSION 

We propose the random field based theory for calcula- 
tion of phase diagram of mixed ferroelectrics and apply it 
to PZT and BZT materials. It has been shown that pro- 



posed theoretical approach describes both quahtatively 
and quantitatively the observed phase diagrams, includ- 
ing relaxor behaviour in BZT at a; > 0.3 (the measured 
fraction x w 0.27). This discrepancy may be related to 
the model assumption that electric dipoles of titanium 
and zirconium ions are the main sources of random fields. 
Zirconium ions can be considered as dilatation centers or 
elastic dipoles which are known to destroy a long-range 
order leading to relaxor properties appearance |j]. The 
calculations of latter property and of the contribution of 
oxygen octahedra tilting to the polarization of PZT are 
in progress now. 
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